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回 热 不 完全 性 以 及 冷 热源 间 存 在 的 热 漏 等 不 可 逆 因 素 ， 继 而 得 到 了 该 模型 的 效率 、 功 率 以 及 生态 学 性 能 函数 。 考 虑 到 三 个 优 
化 指标 不 可 能 同时 达到 各 自 的 最 优 值 ， 故 采用 多 目标 遗传 算法 对 其 进行 三 目标 同时 优化 ， 并 用 TOPSIS 法 进行 决策 并 将 最 终 的 
结果 与 单 目标 优化 进行 了 对 比 。 
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Abstract: This paper presents an investigation of a finite time thermodynamic analysis of a Stirling engine model 
based on the linear phenomenological heat transfer law. Some irreversibilities such as heat resistance, imperfect 
regeneration and heat leak between heat reserviors are considered. As a result, the thermal efficiency, output power 
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a Ahmadi 等 人 [22-2 做 了 大 量 的 工作 ， 用 多 目标 遗传 
HU 局 算法 对 斯 特 林 发 动机 进行 功率 、 效 率 、 热 经 济 学 函 


斯 特 林 发 动机 是 一 种 结构 简单 、 安 全 可 靠 、 可 
以 非常 有 效 地 利用 能 源 和 对 环境 友好 的 热机 ， 并 且 
由 于 斯 特 林 发 动机 可 以 在 温差 相对 很 低 的 冷 热源 间 
做 功 ， 故 其 可 以 燃烧 任何 品质 的 任何 能 源 作 为 其 热 
源 叫 。 因 此 ， 斯 特 林 发 动机 近年 来 引起 了 越 来 越 多 
的 关注 ， 并 被 认为 是 未 来 极 具 潜 力 的 一 种 热机 。 过 程 组 成 ， 如 图 1 所 示 ， 具 体 为 1) 工 质 在 压缩 

斯 特 林 发 动机 在 1816 年 由 罗伯特 。 斯 特 林 发 腔 放 热 压缩 (2) 工 质 进入 回 热 器 进行 等 容 吸 
HH. Blanik 中 等 人 建立 了 内 可 逆 斯 特 林 发 动机 模型 HA; (3) 工 质 进入 膨胀 腔 进行 等 温 膨胀 ， 对 外 做 
并 对 功率 进行 优化 。Thombare ”等 人 Kongtragool Th; (4) 工 质 进入 回 热 器 进行 等 容 放 热 。 斯 特 林 

和 Formosa EAO 征用 等 温 模 型 研究 了 不 完全 回 热 机 的 工 质 不 断 地 进行 这 四 个 过 程 完成 循环 做 功 。 
的 斯 特 林 机 中 死 容积 对 性 能 的 影响 。Tlili AE! 
建立 了 斯 特 林 机 不 可 道 循环 模型 ， 研 究 了 回 热 效率 
和 和 死 容积 的 影响 。 自 从 Curzon and Ahlborn [6 研究 

> 了 存在 有 限 速率 传 热 的 Carnot 热机 ， 开 创 了 有 限 
人 时 间 热 力学 的 先河 之 后 ， 有 限时 间 热 力学 的 方法 也 
CO 开始 被 用 于 分 析 斯 特 林 发 动机 。Ladas ”中 等 人 将 有 
O 限时 间 热 力学 的 方法 应 用 在 斯 特 林 发 动机 上 并 讨论 
© 了 换 热 时 间 对 功率 和 效率 的 影响 。Costea 和 
re Petrescu Al 9 用 基于 有 限 速度 热力 学 第 一 定律 在 实际 情况 下 ， 由 于 热 〈 冷 ) 源 和 工 质 之 间 存 


数 、 炉 产 率 、 压 降 等 进行 多 目标 优化 。 


1 斯 特 林 发 动机 热力 学 分 析 


经 典 的 斯 特 林 循环 由 两 个 等 温 过 程 和 两 个 等 容 


1 斯 特 林 循 环 示意 图 
Fig.1 Stirling cycle diagram 


O 对 斯 特 林 发 动机 进行 了 性 能 分 析 。Senft nq 考虑 了 热 阻 ， 工 质 和 热 ( 冷 ) 源 并 不 是 同一 个 温度 。 此 


在 

有 限 传 热 、 机 械 损 失 和 内 部 热 损 失 的 情况 下 斯 特 林 ”外 ， 实 际 的 回 热 器 并 不 能 有 百分之百 的 回 热效率 ， 
© 机 的 理论 模型 。Chen ”等 人 5 用 有 限时 间 热力 学 的 ”会 有 一 定 的 回 热 损 失 存在 ， 所 以 工 质 在 经 过 回 热 过 
“方法 建立 了 太阳 能 斯 特 林 机 的 理论 模型 ， 得 到 了 系 。 ” 程 之 后 并 不 能 达到 理想 点 。 如 图 2 所 示 。 

之 统 的 效率 上 限 。Wu 等 人 03 用 有 限时 间 热力 学 的 方 

法 分 析 了 在 回 热 不 完全 的 情况 下 斯 特 林 机 的 模型 并 
得 到 了 净 输 出 功 和 效率 的 关系 。Kaushik 等 人 03 1 
用 有 限时 间 热 力学 的 方法 考虑 了 有 限 热 容 热源 的 情 
况 下 考虑 回 热 损 失 的 斯 特 林 模 型 ， 并 分 析 了 各 参数 
对 功率 和 效率 的 影响 。 

研究 人 员 对 斯 特 林 机 做 了 大 量 的 优化 工作 。Li 


2 考虑 实际 情况 的 斯 特 林 机 循环 示意 图 


等 人 05] 用 有 限时 间 热 力学 的 方法 对 太阳 能 斯 特 林 发 Fig.2 Stirling machine cycle diagram of the actual 
动机 的 功率 进行 优化 ， 并 得 到 相应 的 效率 。 situation 
Angulo-Brown 06] 提出 了 一 种 生态 学 函数 作为 热机 在 该 系统 中 ， 一 个 循环 内 热源 对 斯 特 林 发 动机 


的 评价 指标 ， 使 热力 循环 的 功率 与 功率 耗 散 达到 最 内 工 质 的 传 热量 为 
佳 的 折衷 。Yan07 对 Angulo 的 指标 进行 改进 ， 提 出 Ei 

使 用 环境 温度 代替 冷 源 温度 更 为 合理 。Long U8) i ou.-o | 二 } (1) 
究 了 具有 非 等 温 过 程 一 般 化 热机 在 生态 学 准则 下 的 

效率 限 。He 与 Tyagi 等 人 09 20 用 生态 学 函数 作为 ”其 中 a 为 热 端 换 热 器 的 传 热 系数 ， 7 为 热源 温 
目标 函数 对 不 可 逆 斯 特 林 机 进行 优化 。Ust PURE ” 度 ， 有 为 工 质 在 膨胀 腔 的 温度 ， 厂 为 一 个 循环 内 工 


了 新 的 生态 学 函数 一 一 ECOP 来 对 热机 进行 评价 ， 该 
评价 指标 表示 为 功率 和 功率 耗 散 之 比 。 近 年 来 ， 涌 
现 了 一 些 对 斯 特 林 发 动机 进行 多 目标 优化 的 研究 。 


质 在 膨胀 腔 的 时 间 。 同 样 可 得 一 个 循环 内 冷 源 与 斯 
特 林 发 动机 内 工 质 的 换 热 量 为 


1 1 
On = a z 一 =) b (2) 


其 中 6 为 冷 端 换 热 器 的 传 热 系数 ，7 为 冷 源 温 
度 ， 刀 为 工 质 在 压缩 腔 的 温度 ， 乌 为 一 个 循环 内 工 
质 在 压缩 腔 的 时 间 。 

在 3-4 过 程 中 ， 工 质 进行 的 是 等 温 膨 胀 过 程 ， 
所 以 热源 传递 给 工 质 的 热量 等 于 工 质 对 外 做 功 量 ， 
故 有 


Q., =nRT, Ind (3) 
其 中 z 为 工 质 的 摩尔 数 ， 为 理想 气体 常数 ， 4 为 
气体 膨胀 压缩 比 。 在 斯 特 林 循 环 中 ， 有 
ae 
V, OV; 

同样 的 ， 在 1-2 的 过 程 中 ， 工 质 对 冷 源 的 放 热 
量 等 于 外 界 对 工 质 的 做 功 量 


(4) 


Qp =nRT, Ind (5) 
所 以 ， 工 质 在 膨胀 腔 和 压缩 腔 的 时 间 分 别 为 
nRT Ind 
j œ (Ta -T,) 


回 热 过 程 2-3 和 4-1 均 为 等 容 过 程 ， 无 做 功 ， 
所 以 回 热 器 与 工 质 的 换 热 量 等 于 工 质 内 能 的 改变 


2 =c,m(T, =D) (7) 


Q, = c,m(T, T T,) (8) 


其 中 c HTRRWERAZ, mAs, mM 刀 分 

别 为 等 温 过 程 3-4 和 1-2 中 工 质 的 温度 。 

在 实际 过 程 中 ， 由 于 存在 回 热 损 失 ， 工 质 在 回 
过 程 中 达 不 到 温度 A 7/， 而 是 m 和 Ty ， 传 

a JE Q M On, ME Ay Mar» MAA 


2 =cm(T,.—T,) (9) 
Qr =cm(T, —T,) (10) 
假设 2-3 和 4-1 过 程 的 回 热 损 失 系 数 分 别 为 4 
1 和 Ho, 于 是 有 
3 = (1— 4) Oo, (11) 


ar = (1-4) Ou, (12) 


联 立 (3) - (8) 式 可 得 


Ty =(1- 4 )T, + 41, (13) 
T, = 41, +=) (14) 

在 回 热 过 程 中 ， 为 求 回 热 时 间 ， 定 义 回 热 时 间 
系数 
et ae (i=1,2) (15) 
dt 
其 中 ，“+” ” “~” 分 别 对 应 升温 和 降温 过 程 ， 六 1 


对 应 升温 过 程 ， 庆 2 对 应 降温 过 程 。 当 温度 变化 不 
均匀 时 ，k; 代 表 工 质 温 度 的 平均 变化 率 。 继 而 可 得 
质 在 回 热 器 里 的 时 间 分 别 为 


T, -T. 
by = (16) 
ly = 2 = (17) 
k, 
其 中 及， 和 九 ， 是 考虑 了 不 完全 回 热 损 失 的 工 质 在 


回 热 过 程 后 所 达到 的 温度 。 

由 于 回 热 损 失 的 存在 ， 仅 凭 回 热 器 的 回 热量 无 
法 使 得 工 质 回 到 理想 循环 该 有 的 温度 ， 故 需要 靠 热 
源 和 冷 源 额外 的 放 热 ， 所 以 在 一 个 循环 中 ， 热 源 对 
工 质 的 实际 放 热量 和 冷 源 对 工 质 的 实际 吸 热量 为 


Q, = OQ + Qs, =(1- 


4 )c,m(T,-T,)+nRT, nda 


(18) 
O, = Oy + OQ = (1-4 )c,m(T, -T,)+nRT, Ind 
相应 的 时 间 分 别 为 
a,(T,-T,) 
ea eT hee) (21) 


p= 
< œ (T,-T,) 


在 实际 的 斯 特 林 发 动机 运行 中 ， 热 源 和 冷 源 之 
间 有 热 漏 的 存在 ， 在 一 个 循环 中 热 漏 的 总 量 为 


1 1 
Oreak = Cs 2- (22) 


HH, CERIA EA, 
期 ， 根 据 式 子 (16) (17) (20) (21) 可 知 ， 


r 为 循环 周 
r 可 表达 为 


T =b thy thoy thay 
(1-a )em(T, -T,)+nRT, in 2 


a (T -T,) 
Cems )em(T, —T,)+nRT, Ind 
a, (T, -T,) 
-r B-T 
k k, 


(23) 


所 以 ， 由 于 热 漏 的 存在 ， 热 源 实际 的 放 热量 可 
以 表示 为 
Op =Q, + Denk =(1- 4, )e,m(T, -f,) 


snk miia 
L H 


(24) 
同 理 可 知 冷 源 实际 的 吸 热量 可 以 表示 为 


O, = Q, + Qeat = (1-4, )e,m(T, -T,) 


4+nRT, nA +C; +-+] 
在 该 模型 中 ， 斯 特 林 发 动机 的 功率 为 
P=Q,-Q,=nR(T,-T,)Ind (26) 


所 以 效率 可 以 表示 为 
77 = P 
Q, 


=nR(T,-T,)na/[(1-4,)e,m(T, -T,) (27) 


+nRT, Ina +C, Ca 


L H 
由 Ust 提出 的 生态 学 性 能 优化 准则 
ECOP= 2 (28) 
o0 


的 最 大 化 用 来 表示 功率 和 功率 耗 散 的 最 佳 折 袁 。 其 
中 环境 温度 DRR o 的 乘积 ho 被 称 为 “也 
在 此 发 动机 的 模型 中 ， 生 态 学 性 能 函数 可 以 表 


en 


示 为 


=nR(T,-T,)Ina 
C,(T, -T,) A+ 4cm(T, -T,)+nRT, Ind 


(T 

(T T, 

a a sc 
Ty 


2 多 目标 优化 


在 实际 优化 问题 中 ， 经 常 存 在 多 个 优化 指标 ， 
因为 不 可 能 同时 让 所 有 指标 达到 各 自 的 最 优 值 ， 所 
以 往往 出 现 一 个 优化 目标 达到 最 大 值 时 ， 其 它 指标 
并 不 是 最 优 甚至 是 最 劣 值 。 本 文采 用 多 目标 遗传 算 
法 对 斯 特 林 发 动机 的 功率 (26) 、 效 率 (27) 以 及 
生态 学 性 能 准则 〈29) 进行 优化 。 发 动机 工 况 参数 
选择 如 下 : 

工 质 摩尔 数 大 204 mol， 工 质 质量 1816 g， 工 质 等 
容 比热容 c,=3. 214 J/(g。K) ， 理 想 气体 常数 
8.314 J/(mol + R， 回 热 损失 系数 Li=0.3，L4 
2=0.2， 回 热 器 中 工 质 温 度 平 均 变 化 率 f=fp=500 
K/s， 换 热 系数 W/(K。s)， 环 境 温 度 H=300 K, Æ 
源 与 冷 源 间 热 导 率 c=10 WK -e s). 

采取 四 个 决策 变量 p Tr Tr Thy TIRAR 
工 质 在 高 温 膨胀 时 的 温度 、 工 质 在 低温 压缩 时 的 温 
度 、 热 源 温度 、 冷 源 温度 。 决 策 变量 的 约束 条 件 如 
下 所 示 


600 < T, <1000 K (30) 
300 < T, < 600 K (31) 
1000 < T, <1200 K (32) 
280 < T, <300 K (33) 


在 多 目标 优化 问题 中 ， 得 到 的 帕 累 托 最 优 解 是 
一 个 解 集 〈 帕 累 托 前 沿 ) 。 在 得 到 帕 累 托 最 优 解 集 
之 后 ， 需 要 引入 决策 方法 从 中 选取 最 终 的 唯一 结 
果 。 本 文 将 采用 TOPSIS 法 对 帕 累 托 最 优 解 集 进 行 
决策 。 

在 进行 决策 前 ， 需 要 对 不 同 量 纲 的 目标 i 


`N 
p| 


量 纲 化 。TOPSIS 法 使 用 欧式 无 量 纲 化 的 方法 。 

欧式 无 量 纲 化 的 方法 中 ， 帕 累 托 前 沿 上 各 个 不 
同 的 点 组 成 的 目标 矩阵 可 以 表示 为 f;; ， 其 中 了 代 
表 位 于 帕 累 托 前 沿 上 的 每 一 个 点 ，J 代表 目标 空间 
上 的 每 个 点 ， 因 此 可 以 将 无 量 纲目 标 fj 定义 为 : 


F, 
n y 
i ; 


pel) 


TOPSIS 决策 方法 中 ， 定 义 了 理想 点 和 非 理 想 
点 。 理 想 点 代表 着 该 点 上 的 每 个 目标 函数 都 是 最 优 
的 ， 而 在 实际 的 多 目标 优化 问题 中 ， 不 可 能 存在 这 
样 的 理想 点 。 同 样 的 ， 非 理 想 点 代表 着 该 点 上 的 每 
个 目标 函数 都 是 最 劣 的 ， 实 际 问 题 中 也 不 存在 这 样 
的 非 理 想 点 。 因 此 ， 理 想 点 和 非 理 想 点 并 不 在 帕 累 
托 前 沿 上 。 在 TOPSIS 决策 方法 中 ， 对 于 经 过 欧式 
无 量 纲 化 的 目标 ， 每 一 个 帕 累 托 前 沿 上 的 点 与 理想 
点 和 非 理想 点 的 距离 被 定义 为 : 


n 2 
d= Ère) (35) 


j=l 


(34) 


R 2 
d, = [Stn — ppor) (36) 


Jal 


dj, 代表 帕 累 托 前 沿 上 的 点 与 理想 点 的 距离 ，g;。 代 
表 Pareto 前 沿 上 的 点 与 非 理想 点 的 距离 ， 代表 目 
PRET BL, Fidel 和 Fy0r1oe971 分 别 代表 第 7 个 目标 


ChinaXiv 合 作 期 刊 


在 单 目 标 优化 条 件 下 的 最 优 值 与 最 劣 值 。 
定义 参数 Cl; 
d, 


Cl, = 一 -一 37 
d,, +d. = 
而 拥有 最 小 值 C 的 解 被 选 作 最 终 解 ， 因 此 最 终 解 

的 指数 7 可 有 如 下 表示 
inna =i € min (CL );i =1,2,.. , m (38) 


3 结果 与 讨论 
功率 (26) 、 效 率 OD 和 生态 学 性 能 指标 
(29) 在 本 优化 中 被 同时 最 大 化 ， 基 于 多 目标 遗传 
算法 得 到 的 帕 累 托 最 优 解 集 如 图 3 所 示 。 图 中 的 点 
集 即 为 所 求 的 帕 累 托 前 沿 。 经 过 拟 合 之 后 ， 可 以 得 

到 


ECOP =1.16x107 


2 2 
_1 f P+2.61x10° | 1 {743.6110 
2.60x10° J 2 \ 9.43x107! 


-1.17x102xe ° (39) 
拟 合 优 度 为 : 
R? = 0.99867 (40) 


将 TOPSIS 决策 得 到 的 结果 与 单 目标 遗传 算法 
优化 的 结果 进行 对 比 ， 如 表 1 和 表 2 所 示 。 


0.45 


P (CW) 


3 多 目标 遗传 算法 优化 结 
Fig.3 Optimization Results of Multi-objective Genetic Algorithm 


表 1 多 目标 优化 与 单 目 标 优化 结果 对 比 1 
Table 1 The comparison between results of multi-objective optimization with single objective optimization 
目标 函数 变量 

P (kW) | N(%) | EOP | 7,(K) TK) 7; (K) 7, (K) 
单 目标 优化 功率 48. 39 32. 03 0. 67 865. 2 482.9 1200. 0 280. 0 
优化 优化 效率 29. 52 40. 01 1. 04 988. 2 341.8 1169. 7 280. 0 
优化 生态 学 20. 39 37. 73 1.15 910.8 351.7 1017. 3 300. 0 
多 目标 优 (TOPSIS) 42. 37 37. 60 0. 94 948. 8 414. 3 1198.9 288. 1 


表 2 多 目标 优化 与 单 目 标 优 化 结果 对 比 2 


Table 2 The comparison between results of multi-objective optimization with single objective optimization 2 


VS 功率 优 -12.44% +17.39% +40.30% 
VS 效率 优 +43.53% -6.02% -9.62% 
VS 生态 学 函数 优化 | +107.80% -0.34% -18.26% 
由 对 比 表 1 MK 2 可 以 看 到 ， 虽 然 单 目标 优化 的 结果 虽然 不 是 各 项 指标 中 最 大 的 ， 但 是 却 能 达到 
的 结果 可 以 使 得 该 单 目 标 取得 最 大 值 ， 但 是 相应 个 综合 最 优 的 结果 。 例 如 ， 相 比 于 单 目标 优化 功 
的 ， 其 它 指标 却 要 低 得 多 ， 例 如 单 目标 优化 功率 的 率 ， 虽 然 少 了 约 6 kW 的 做 功 量 ， 但 是 却 有 着 约 


时 候 ， 虽 然 功率 是 最 大 的 ， 但 是 效率 和 生态 学 性 能 5.6 % 的 效率 提升 和 约 0. 3 的 ECOP 的 提升 。 相 比 于 
指标 都 是 最 低 的 。 由 表 1 还 可 以 看 到 ， 多 目标 优化 单 目标 优化 效率 ， 效 率 昌 然 低 了 2.4 %， 但 是 功率 


I 


却 有 着 12. 85 kW 的 提升 。 同 样 的 ， 相 对 于 单 目标 优 
化 ECOP 指标 低 了 0. 2， 但 是 多 目标 优化 的 结果 却 有 
着 22 kW 的 提升 。 优 化 的 结果 显示 ， 采 用 多 目标 优 
化 可 以 对 斯 特 林 发 动机 进行 效率 、 功 率 以 及 生态 学 
性 能 指标 有 着 更 合理 的 优化 效果 。 


4 结论 

本 文 对 斯 特 林 发 动机 进行 有 限时 间 热 力学 的 分 
析 ， 考 虑 了 热源 与 工 质 传 热 的 不 可 逆 性 、 回 热 损 失 
和 热 漏 的 影响 ， 再 对 该 模型 求 出 其 功率 、 效 率 以 及 
生态 学 性 能 指标 。 考 虑 到 单 目 标 优化 的 结果 会 使 其 
它 指标 付出 很 大 的 代价 ， 不 具有 实际 指导 意义 ， 故 
采用 遗传 算法 对 这 三 个 指标 进行 多 目标 优化 。 在 多 
目标 遗传 算法 进行 优化 得 到 帕 累 托 最 优 解 集 之 后 ， 
采用 了 TOPSIS 决策 方法 进行 决策 ， 得 到 的 结果 与 
单 目标 优化 的 结果 进行 对 比 。 从 对 比 的 结果 可 以 看 
到 ， 采 用 多 目标 遗传 算法 对 斯 特 林 发 动机 进行 优化 
的 结果 更 具有 实际 指导 意义 。 
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